Finite scale singularity in the renormalization group flow of a reaction-diffusion system 
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We study the nonequilibrium critical behavior of the pair contact process with diffusion (PCPD) 
by means of nonperturbative functional renormalization group techniques. We show that usual 
perturbation theory fails because the effective potential develops a nonanalyticity at a finite length 
scale: Perturbatively forbidden terms are dynamically generated and the flow can be continued once 
they are taken into account. Our results suggest that the critical behavior of PCPD can be either 
in the directed percolation or in a new (conjugated) universality class. 



Reaction-diffusion systems involving one particle 
species (also known as branching and annihilating ran- 
dom walks (BARW)) are stochastic out of equilibrium 
systems important from both a phenomenological and a 
theoretical viewpoint. They consist of identical particles 
A, diffusing on a d-dimensional lattice, that can branch 
{nA —> (n + p)A) or annihilate (nA —> (n — q)A). The 
competition between these two types of reaction is gener- 
ically responsible for the existence of transitions between 
an active phase where the density of particles is finite, 
and an "absorbing" phase where all particles, and thus 
all fluctuations, have disappeared. Such models provide 
the building blocks of a large variety of applications and 
models in physics and beyond, and are therefore of fun- 
damental importance pQ. They also have the advantage 
of providing a relatively simple theoretical framework for 
the study of the different universality classes of absorbing 
phase transitions. 

Our understanding of out-of-equilibrium critical phe- 
nomena in general and absorbing phase transitions in 
particular has benefited from perturbative approaches 
O [3] , but important advances were brought recently by 
the application of nonperturbative renormalization group 
(NPRG) methods [2HSH]- For the two prominent cases of 
BARW, A -> 2A, 2A -> which represents the directed 
percolation (DP) class, and A — > 3A, 2A — > which be- 
longs to the parity-conserving, or generalized voter class, 
the success of the NPRG owed to the presence of non- 
perturbative features. 

The case of the "pair contact process with diffusion" 
(PCPD) has largely resisted analysis so far [8] but was 
not studied with NPRG methods. The simplest BARW 
model in this class consists of reactions 2A — > 3A and 
2A —> with rates a and A. (A limiting reaction such 
as 3^4 — > 0, with rate A', is actually needed to ensure a 
finite density active phase [22].) The distinctive feature 
of the PCPD is that two particles must meet to trig- 
ger branching. On general grounds, this is not expected 
to be a relevant ingredient defining universality classes, 
hence the interest raised by results obtained so far on the 
critical behavior of PCPD: It has been intensively stud- 



ied numerically in d = 1 but remains unclear because of 
the presence of slow dynamics and/or strong corrections 
to scaling 13J . The debate, ongoing still recently, is to 
know whether PCPD belongs to the DP universality class 
[SHU] or not [SHU]. 

Even the status of d c , the upper critical dimension of 
PCPD, is unclear: numerically, d = 3 seems beyond it 
|20) . but in d = 2, the presence of large corrections to 
scaling is difficult to disentangle from logarithmic terms 
preventing clear conclusions in spite of indications of 
mean- field behavior [51] [22] . In perturbation theory the 
RG flow of PCPD goes to the Gaussian fixed point for 
d > 2 and sufficiently small coupling constants, whereas 
it blows up at a finite scale for larger couplings or for 
d < 2 9J. This suggests d c = 2. Note that the ex- 
plosive flow forbids the exploration of the long-distance 
physics of the model. This is also known to occur in 
quantum chromodynamics (at the confinement scale), in 
the O(N) nonlinear sigma model (at the scale of the cor- 
relation length) |10) . and in pinned elastic manifolds (at 
the Larkin length) [TT]. In this last case, NPRG methods 
at the functional level allowed to treat the problem [T2"] . 

In this Letter, we examine the PCPD field theory in 
the light of the NPRG, explain why perturbation theory 
fails, and how to avoid its problems. We show that the 
potential in the running effective action develops a sin- 
gularity at a finite scale, signalling that couplings that 
are perturbatively forbidden are dynamically generated. 
Once taken into account the RG flow can be continued 
and a fixed point can be found. Our results suggest that 
the critical behavior of the model is either in the DP class 
or possibly in a new class characterized by a "conjugated" 
symmetry, and that d c — 2 only at small coupling. Our 
study indicates that NPRG is a powerful tool for dealing 
with similar situations beyond reaction-diffusion systems. 

The field theory associated with PCPD. By using the 
usual Doi-Peliti formalism [22] it is possible to derive the 
action associated with PCPD from first principles: 

S= f [0(9^-i?V 2 0)+0 2 ( ffl + . 92 2 + ff3 3 ) 
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where 4> and <j> are complex conjugates, x = (t, x), J = 
J d d x dt, and g\ = 2A — a, .92 = A — 2a, g 3 = — a. Within 
perturbation theory, one finds that initializing the RG 
flow with the (bare) couplings of the action S, Eq.([T]), 
the running coupling g 2 (k) (k being a momentum scale) 
diverges at a finite scale k c for d < 2 which bars from 
exploring scales below k c . Moreover, the {<j>(f)) response 
function does not receive any loop correction and thus the 
dynamical exponent z remains equal to 2 which is clearly 
invalidated by numerical results. For d > 2, all couplings 
are irrelevant around the Gaussian fixed point which is 
thus (locally) attractive. (Hence the conclusion that d c = 
2.) We now present a calculation where the flow no longer 
diverges but develops a singularity that can naturally be 
taken into account at the price of working functionally. 

The nonperturbative renormalization group. NPRG 
follows Wilson's idea [23] of partial integration over fluc- 
tuations. It builds a one-parameter family of models in- 
dexed by a momentum scale k such that fluctuations are 
smoothly included as k is lowered from the inverse lat- 
tice spacing A down to k = where they have all been 
summed over. To this aim, we add to the original action 
a momentum-dependent mass-like term |25l I26j 



(2) 



where q = (w, q), q = \q\, i = 1, 2, <j> x = 4>, (f> 2 = (j) and 
repeated indices are summed over. With, e.g., [R k ]\ 2 — 
[Rkhi = (k 2 - q 2 )9(k 2 - q 2 ) and [R k ] xl = [R k } 22 = 0, 
the fluctuation modes <fri(q > k) are unaffected by AS k 
while the others with q < k are essentially frozen. The 
fc-dependent generating functional of correlation and re- 
sponse functions thus reads 



Z k [Ji,J 2 ] 



D<t>D<t)e 
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The effective action r^^], where ipi = ip = (4>i), 1P2 = 
ip = (4> 2 ), is given by the Legendre transform of W k = 
logZfe (up to the term proportional to R k ): 

T k [ip,^}+W k = f Jrfi- f i? fe (gh%)^(-q). (4) 

The two-point functions can be computed from T k by 
differentiating: 



s 2 r k 



(5) 



5Vu(Xl)^is ( x 2) 

and the exact flow equation for T k [i/j, ^] reads [25] : 

d k T k = ^Tvj d k R k . G k with G k = [T {2) + R k ]-\ (6) 

When k decreases from A to 0, T k varies between the 
(bare) action: T k= \ — S, and the full effective action: 



r k=0 = T. Solving the flow equation ^ is thus equiv- 
alent to solving the model. This is however impossible 
to do exactly and approximations must be made. We 
perform here the local potential approximation (LPA) 
which is known to work well for the determination of the 
critical behavior of models either at equilibrium or out- 
of-equilibrium [3 H2 HZ] : 



r . t^LPA 
1 k -> 1 k 



(^(d t -DW 2 )tP + U k (^,i>)). (7) 



Substituting Eq.([7]) into Eq.(|6|, choosing the R k function 
described above and integrating over q, we find : 



d k U k = C d 



k 2 + U, 



(1,1) 



- 1 



(0,2)^(2,0) 



(8) 



where C d = 4k 2+d [2 d+1 ir d / 2 dT(d/2)}- 1 and the upper 
indices code for derivatives in ip and ip. It is easy to ver- 
ify that performing a field expansion of Uki^),^) around 
(0, 0) leads, for the couplings in front of the monomials 
ip m ip n , to the same difficulties as encountered perturba- 
tively: the flow blows up at a finite scale for d < 2. We 
thus need to work functionally. However, we can further 
simplify our ansatz Eq.([7| by performing an expansion 
of U k {if),ip) in if) only while remaining functional in the 
■0-direction. We thus replace U k (ip,'ip) by 



JV 



(9) 



and we typically truncate the sum at order N = 3 or 
4. Our approximation scheme is therefore based on the 
assumption that the LPA is sufficient, and that the non- 
trivial features of the model lie in the ip direction (which 
is not spoiled by the expansion in ?/>). Inserting Eq.|9| 
into Eq.|8]) we find the flow equations for the functions 
V n k at order N = 3 : 



Vi = BV 2 V{' 
'1 



v 2 
v 3 



B 



B 



V" (AV 2 (3AV2V" - 8V£) +6V3) + V 2 V 2 " 

l - (A (8V 2 ' (3V[' (AV 2 V 2 ' - V 3 ) - V 2 V 2 ") 
+V 2 V(' (A (V(' (AV 2 (5AV 2 V" - 2AV 2 ') + 18V 3 ) 



(10) 
(11) 



6V 2 V 2 ") - 12F 3 ')) + 6V 3 V 2 ") + V 2 V 3 " 



(12) 



where, for simplicity, we have omitted the index k in the 
functions V n , kl and V£ are the derivatives of these 
functions with respect to if>, V n = kd k V n , A — (k 2 +V{)~ 1 
and B = (k d+2 Av d /d)A 3 \A~ 1 \. The initial conditions of 
these flow equations are provided by the action S with 
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FIG. 1: (color online) Left: Flows in d = 1 of the linear 
couplings a n (k) = V^ k (0) in the semi-functional approach, 
Eq.(|9|. Right: Same when Uk is expanded in both ip and ip 
and infinitesimal linear couplings are considered (<ri(fc = A) = 
0, (72 (k — A) — 10~ 8 , all other rates are of order 1). Only 
terms at most quartic in the fields have been retained. The 
RG flow of all couplings remains finite but shows an abrupt 
change around k c ~ 0.6. 



4>i —> ipi- They impose that Vn: V£ k= ^{ip = 0) = 0: the 
bare potential has no linear term in <p. Perturbatively, 
it is easy to show that these terms cannot be generated 
in Uk since all Feynman diagrams involve at least two 
incoming particles and thus two fields ip. 

We have numerically integrated the coupled flow equa- 
tions of the functions V n ,k{i>) with N 
for the N 



4 (see Eqs.([10|) 
3 case) together with the initial conditions 
provided by S, Eq.|l]). In the early stage of this flow, the 
linear term of each of these functions remains identically 
zero as naively expected. However, at a finite scale k c 
which is typically the scale where the perturbative flow 
blows up, a linear term is generated in all these functions 
(Fig{T]): the potential Uk develops a corner at ip = and 
its analytic structure is changed below fc c . 

A detailed study of the emergence of the linear terms 
reveals that for k > fc c , a boundary layer appears in 
the V nj k{ip) functions such that in the inner part of the 
layer, that is, at small ip, these functions are expandable 
around ip = (for ip > 0) and start quadratically in ip. 
In the outer region, a linear part appears in the K l ,fe(^/') 
functions. As k approaches k c from above, the width of 
the layer decreases and vanishes at k c , leaving the linear 
term as the dominant term around ip = (FigM. Below 
fe c , the linear terms remain present in the vicinity of ip = 
and the flow can be continued all the way to k — 0. We 
have checked that this emerging scenario holds at order 
N = 3 and 4 of the i/5-expansion. Converged, higher- 
order results are unfortunately difficult to obtain 33 . 

We have confirmed our scenario by considering the 
usual set of reactions 2 A -t 3 A, 2 A -t and 3 A ->• 
complemented by A — > 2A and A — > with infinitesimal 
rates. Expanding the potential Uk in both ip and ip we 
now have (at least) two new couplings that are linear in 
ip: <7i(k)ipip and o~2{k)ipip 2 with a\^,(k — A) infinites- 
imal. If the limit a\^(k = A) — > were regular then 
0i,2 (&) would start playing a significant role in the RG 
flows of the other couplings only below a very small scale 
k that would go to as a\^{k = A) — > 0. We find on 
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FIG. 2: (color online) Boundary layer for V^fc (?/>)> whose be- 
havior around the origin is displayed for k ~ k c and several 
values of k above k c . For all k > k c , V% t k starts quadratically 
at ip = whereas a linear term appears at k = k c . Inset: 
width of the boundary layer as a function of logfe/A. 



the contrary that while for k just below A the flows of all 
the couplings are indeed almost insensitive to o~\£ (when 
they are initially extremely small), this is no longer the 
case for k ~ k c since the dramatic increase of g2(k) makes 
0i,a(fc) grow abruptly around k c , independently of their 
initial smallness, see Fig. [T] For k close to k c , the back 
reactions of cri^fe) on the flows of all the other couplings 
start to be significant and eventually modify them com- 
pletely since o~\p are the most relevant couplings. The 
RG flow is no longer singular (32 (k c ) remains finite) but 
lives below k c in a larger functional space involving the 
couplings linear in ip. This result is fully consistent with 
what is found in the functional viewpoint. 

Criticality in d = 1. Our two different approaches 
both conclude that terms linear in ip are generated be- 
low a nonuniversal scale k c in d = 1. We can therefore 
consider the field theory obtained just below k c as a new 
field theory that can be studied per se, the difficulty be- 
ing that its action is non polynomial since the functions 
V n ,k{ip) are not. In a perturbative analysis only the terms 
of lowest degrees in ip and ip would be retained in the bare 
action, i.e. ipip, ip 2 ip and ipip 2 . Depending on the relative 
sign of the two cubic terms, this action, truncated at or- 
der three, exhibits one of the following symmetries (after 
a trivial rescaling of the fields): ip(t) ^ ±ip(—t). The 
minus sign corresponds to the cubic terms having oppo- 
site signs. This is the "rapidity" symmetry defining the 
DP class. The other sign defines a new, "conjugated", 
symmetry. We call the corresponding class DP'. It is 
easy to show in our framework, Eq.([8|, that when only 
the above terms are kept in Uk, only two nontrivial fixed 
points exist and that they show either one or the other 
of the above symmetries (we call them DP and DP'). 
No such result exists beyond this simple truncation. To 
the best of our knowledge, the DP' symmetry has never 
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FIG. 3: (color online) Left: Flow of the symmetry-testing 
ratio rfe computed in the semi-functional approach with TV = 
4. For k < k c , rk flows toward the DP value r k DP ^ = —1. 
Right: line separating the basin of attraction of the Gaussian 
fixed point and the domain where the singularity is present 
in the (d, (72(A)) plane when Pi(A) = 0. 



appeared nor been studied before. 

Two difficulties appear when studying the existence 
and the nature of the fixed point. First, when using a 
complete field expansion of Uk, it is difficult to initialize 
the flow below k c with what has been found function- 
ally using Eq.(|9| just after the singularity because this 
amounts to projecting functions of tp onto polynomials, 
setting infinitely many terms of high degree to 0. Sec- 
ond, if we work directly with our semi-functional approx- 
imation, Eq.|9|, the roles of ip an d V> are dissymmetric 
which hinders the search of a fixed point exhibiting a 
symmetry that exchanges ip and ip as in DP or DP'. We 
have nevertheless studied the existence of a fixed point 
using both approaches. Within richer and richer poly- 
nomial approximations of Uk (where ip and ip play sym- 
metric roles) initialized with couplings found with the 
semi-functional approximation just after the singularity, 
we have found a range of initial reaction rates leading to 
the DP fixed point. In the other, semi- functional approx- 
imation truncated at order N = 4, we have recorded the 



flow of r k = (U { k 2 ' 1] /U { k h2) 



U, 



(1,3) /7T (3,1) 



/U k (calculated at 



ip = ip = 0). This quantity is an indicator of the na- 
ture of the fixed point because if r*. = — 1 the expansion 
of Uk up to order four in the fields exhibits the rapidity 
symmetry (up to a rescaling of the fields). If = 4-1 
the same holds true for the DP' symmetry. Notice that 
just below k c , r k is neither 1 nor —1. We show in Figj3] 
the evolution of computed from Eq.(|9| with N = 4. 
Although the instabilities of our numerical code prevent 
us from finding a true fixed point value for there is 
little doubt that it indeed reaches a plateau at the value 
— 1 so that the expected fixed point should be that of DP 
in agreement with what is found in the polynomial ap- 
proximation. Notice however that if we work with a field 
truncation of Uk, starting at k = A with a very small 
coupling era (fc = A) as explained above, we find, at least 
in the simplest truncation of Uk, either the DP or DP' 
fixed points depending on the initial rates (we have not 
been able to confirm this result with truncations involv- 
ing higher powers of the fields because of the extreme 



sensitivity of the flow at small k on the choice of the 
initial rates). Interestingly, we find that the vicinity of 
both fixed points is reached after a very long RG "time" 
(log fc/A ~ —9) which means that the scaling regime only 
appears at very large lengthscales. This could explain 
why it is so difficult to observe the asymptotic regime in 
numerical simulations. Our general conclusion is there- 
fore that the critical behavior of PCPD in d = 1 should 
be either in the DP or DP' universality class and our best 
results are in favor of DP. 

The upper critical dimension. The determination of d c 
is the same in perturbation theory and in our scheme for 
small initial couplings: for d > 2, the conditions gx,A = 
and <?2,A small make the system critical and the flow is 
driven towards the Gaussian fixed point. However, for 
large enough g 2 , A the flow always becomes singular, even 
when <7i a = 0. We show in Figj3]the basin of attraction 
of the Gaussian fixed point in this case as a function 
of 92, A obtained within perturbation theory. The result 
is similar to what is found in the Kardar-Parisi-Zhang 
equation [2"8] [51]. Note that within perturbation theory, 
criticality is reached when <7i.a = 0. It is no longer clear 
in our approach that this condition is necessary for large 
initial 52, A and d > 2 because the finite scale singularity 
emerges for generic gi.A- Our results suggest that there 
could exist a nontrivial critical behavior at large initial 
92, A even above one dimension that could be in the DP 
(or possibly DP') universality class. 

To conclude, we believe our results represent a break- 
through for the understanding of the critical behavior 
of PCPD. Even though we do not have yet a com- 
plete solution of the problem, we have unlocked an 
heretofore blocked situation and offered new lines of fur- 
ther research. The recourse to functional nonperturba- 
tive renormalization was essential because it allows us 
to address questions such as the generation of linear 
terms. Going beyond our semi-functional approxima- 
tion remains the main challenge that should lead, once 
controlled numerically, to satisfactory results. From a 
physical point of view, a final answer to the question of 
the phase diagram and its accessible fixed points in all 
dimensions remains of course the main goal. But under- 
standing the meaning of the scale fc c is also challenging. 
It could be related to the existence of relevant "elemen- 
tary excitations" made of pairs of particles that would 
involve an intrinsic scale. Disentangling the roles of the 
particles and of the pairs has already been studied by ef- 
fectively taking them into account through the introduc- 
tion of another species B and the reactions: 2A — > B, 
B -> A, B -> 2B, 2B -> B, B -> 0. Such two-species 
PCPD models, which are believed to exhibit the same 
critical behavior as the one species case studied here [31] , 
could be the starting point of more complicated NPRG 
approaches (involving four fields, but possibly deprived 
of singularities in the flow). These endeavors are left for 
future work. 
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